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It has been shown that entanglement distillation of Gaussian entangled states by means of local 
photon subtraction can be improved by local Gaussian transformations. Here we show that a 
similar effect can be expected for the distillation of an asymmetric Gaussian entangled state that 
is produced by a single squeezed beam. We show that for low initial entanglement, our largely 
simplified protocol generates more entanglement than previous proposed protocols. Furthermore, 
we show that the distillation scheme also works efficiently on decohered entangled states as well as 
with a practical photon subtraction setup. 
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Continuous variable (CV) entanglement is a valuable 
resource for many quantum informational protocols pQ- 
|3]. However, the performance of these protocols is often 
limited due to the difficulty in generating CV states with 
a high degree of entanglement. Moreover, even if a large 
degree of entanglement can be produced using a highly 
efficient nonlinear parametric process, the distribution of 
it (e.g. among two parties in a network) will inevitably 
lead to dissipation, rendering the state weakly entangled. 
To improve the performance of quantum information pro- 
cessing, it is therefore important to devise a protocol that 
increases the amount of entanglement between two dis- 
tant parties by means of local quantum transformations 
and classical communication. This can be done by the 
process of entanglement distillation. 

Distillation of non-Gaussian entanglement (either a 
pure or a mixed non-Gaussian entangled state) can be im- 
plemented using simple linear optics [1] as demonstrated 
in [21 inj- On the contrary, the distillation of Gaussian 
states (either pure or mixed Gaussian states) is challeng- 
ing as it will inevitably require some non-Gaussian trans- 
formations [THS] enabled by a very strong Kerr nonlinear- 
ity [l^E], using a non-Gaussian measurement (12^14]. 
using non-Gaussian resources |15| or using a combination 
of photon addition and subtraction fl6'. An intriguing 
scheme for entanglement distillation of Gaussian states 
was suggested by Opatrny et al. [12] and involves local 
photon subtraction of a two-mode entangled state. The 
scheme was recently implemented by Takahashi et al |17| 
(entanglement distillation by non-local photon subtrac- 
tion has also been demonstrated This photon sub- 
traction scheme can however be improved by local Gaus- 
sian operations: Zhang and van Loock [19] showed that 
local squeezing operations prior to photon subtraction 
improve the performance of distillation in terms of pro- 
ducing states with a higher degree of entanglement and 
with higher success rate. Secondly, it was realized by 
Fiurasek [20^ that by using the much simpler operations 
of local phase space displacements, a similar improve- 
ment can be achieved. 



In all these previous proposals on entanglement distil- 
lation, the considered Gaussian entangled state was pro- 
duced by interfering two single mode squeezed states on a 
beam splitter. However, in practice it is much simpler to 
generate entanglement from a single squeezed mode that 
is split on a balanced beam splitter, as was done in [17j . 
In this letter we investigate the displacement enhanced 
entanglement distillation of the simplified Gaussian en- 
tangled state which will be much closer to an experimen- 
tal realization. Surprisingly, we find that in some cases 
the usage of a single mode squeezed beam for the gen- 
eration of Gaussian entangled states, the degree of en- 
tanglement after distillation is higher than if a two-mode 
squeezed state was used. 

We consider the entanglement distillation setup shown 
in Fig. [l] The entangled state is simply produced by 
dividing a single mode squeezed state on a balanced beam 
splitter, and the resulting modes are sent through lossy 
channels to the two sites, denoted A and B. At these 
two sites a Gaussian displacement transformation as well 
as single photon subtractions are applied to distill the 
quantum state. In the following we first consider the 
distillation protocol when the channels are loss-free and 
the photon subtraction is ideal, and secondly we consider 
the more realistic scenario where the channels are lossy 
and the photon subtraction process is non-ideal. 

Our starting point is a single mode squeezed state 
mixed with vacuum on a balanced beam splitter (written 
in the Fock state representation): 

ICi) - f/..(l-7?)^/^f:^(|)"|2.)|0X,l) 

where Ubs is the beam splitter unitary transformation, 
7i — tanh(si) and si is the squeezing parameter. In the 
weak squeezing limit (7^1) the state can be truncated 
to 

ICi) « |00) + ^Q|20) + -^|ll) + i|02)).(2) 
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FIG. 1: a) Schematic of the proposed setup. S denotes the 
single mode squeezed state and D{a), D{l3) are displacement 
operations, b) Two equivalent setups for displacement con- 
trolled photon subtraction where /3' — ^/(l — T)j3 and T is 
the transmission of the beam splitter. 



The entanglement of this state can be quantified by the 
logarithmic negativity (LN) which is defined as the bi- 
nary logarithm of the trace norm of the partially trans- 
posed density matrix (|^i) (Ci|)"^'* ED- The logarithmic 
negativity of the state as a function of the average pho- 
ton number {n) — 7i / (I — 7i ) is plotted in Fig. [2^ by the 
black dotted line. 

The idea is now to increase the entanglement of the 
state, and this can be done by removing the vacuum 
contribution and equalizing the weight of the low order 
excitation terms. For weak squeezing (as in ([2])), it is 
straightforwardly realized that such an equalization can 
be obtained simply by performing single photon subtrac- 
tion on just one of the modes A or B. For single photon 
subtraction of mode A, enabled by the bosonic annihila- 
tion operator aA^ the result is 

ICi) ^ aA®lB\Ci) (3) 

= y(|io> + |oi)), (4) 

which is maximally entangled (in the two-dimensional 
subspace, thus neglecting higher order terms). The re- 
sulting entanglement is plotted in Fig. [2^ with the black 
dash-dotted line, and it is shown that for very low squeez- 
ing degrees the LN reaches the maximum value of 1 for 
a 2D Hilbert space. However, the state should be de- 
scribed in a larger Hilbert space in which the state is 
not maximally entangled. We next consider the simul- 
taneous subtraction of single photons at both sites as 
described by the operator a^® o-b and denoted 2PS (2- 
Photon-Subtraction). The resulting LN is plotted in Fig. 
[2^ by the black dashed line. In contrast to the previous 



protocol (with a single photon being subtracted, IPS), 
this protocol is not very effective for low average photon 
numbers but for higher numbers (larger than about 0.21, 
corresponding to 3.9 dB of squeezing) it becomes more 
effective. 

The entanglement can however be further enhanced 
for all degrees of initial squeezing by applying a Gaus- 
sian transformation prior to photon subtraction. The 
transformation that leads to this enlargement is the sim- 
ple phase space displacements, D{a) = exp(aa^ — a*a^) 
and D(/3) = exp(/3a^ - (3* as), where a and (3 are the 
complex excitations of the displacements. By implement- 
ing these displacements in mode A and B prior to and 
after the two photon subtractions (see Fig.[T]3), the state 
reads 

m = (o/j+|)ioo> + ^o+o)MiJH> 

It is clear from this expression that the vacuum and 
single-photon contributions, which prevent the state from 
being strongly entangled, can be removed by setting 

a = = , (5) 

leading to the state 

l*> = -^Q|20) + i=|ll) + ^|02))+O(7^),(6) 

that contains an entanglement of = 1.54 for low ini- 
tial squeezing degrees. This is not maximally entangled 
in the 3-dimensional Hilbert space due to the unequal 
weights but it is close to - the maximally entangled state 
would have i?jv — 1-585. 

The displacements in ^ are maximizing the LN only 
for low degrees of initial squeezing. For larger squeezing 
degrees, the optimum displacement deviates from the one 
in ([5]) , and thus it is not optimized by the removal of the 
vacuum term. The single-photon components are how- 
ever always cancelled by the choice of a = — /3. The opti- 
mized displacement is plotted in Fig.[2]3 and the resulting 
optimized LN is plotted in Fig. [2^ (both by black solid 
curves), and we clearly see that the pre-Gaussian pro- 
cessing improves the entanglement for all average photon 
numbers [22'. 

For comparison, we also briefly consider the distillation 
of a two-mode squeezed state as was treated in Ref . [20] . 
Here the starting point is 



\C2) AB = I^UHs (7) 

n=0 

where 72 = tanh(s2) and S2 is the two-mode squeez- 
ing parameter. This state, which has an average photon 
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FIG. 2: a) Logarithmic negativity as a function of the average 
number of photons of the initial squeezed state. 1(2)MSV; 
Single (two) mode squeezed vacuum states used initially. 
1(2)PS: Single (two) photon subtraction. D2PS: Displace- 
ment based two-photon subtraction, b) The size of the dis- 
placements (a = — /3) that optimise En as a function of the 
average photon number of the resource states. 



number (n) — 27|/(1 — 7I), can be produced by inter- 
fering and phase locking two single mode squeezed vacua 
on a balanced beam splitter. The implementation of a 
local single photon subtraction, either on one or both 
sites, transforms the entangled state into another entan- 
gled state with higher LN as illustrated in Fig. by 
yellow curves. Note that for identical initial (n), 72 is 
smaller than 71 . 

By displacing the state with excitations of 
aA = — ctB = \/y2i the distilled state reads 
1^72) = 72^^(|10) - |01)) + 0(7!) for weak initial 
entanglement (72 <C 1). This state is also maximally 
entangled in the two-dimensional sub space for low 
initial squeezing as was the case for the IPS single mode 
squeezed state in 

In particular, we note that the degree of entan- 
glement after distilling the two-mode state using the 
displacement-enhanced protocol is lower than the one 
obtained for the single mode scheme if the average pho- 




FIG. 3: Distribution of the different eigenstates of the entan- 
gled state that is produced from a single squeezed mode. The 
distributions are for the initial state (before distillation), the 
single photon subtracted state (IPS), the two photon sub- 
tracted state (2PS) and the displacement-enhanced protocol. 



ton number is low. That is, using the largely simpli- 
fied protocol with a single mode squeezer split on a 
beam splitter (rather than interfering and phase lock- 
ing two single-mode squeezers) the distillation transfor- 
mation produces an entangled state that contains more 
entanglement. However, this conclusion only holds for a 
low average photon number. For higher photon numbers, 
the usage of a two-mode squeezed state results in a state 
with larger entanglement than if a single-mode squeezed 
state was used. But the low-gain regime is often of in- 
terest, since entanglement in this region is fairly easy to 
prepare, manipulate and maintain. 

It is known that the logarithmic negativity is optimised 
when the dominant eigenstates of a given state have equal 
weights as is the case in (|4| and partially in (|6|. To get 
some physical insight into the formation of stronger en- 
tanglement, in Fig. [3]we plot the weights of the different 
eigenstates resulting from the different distillation pro- 
tocols using a single mode squeezed state. Due to the 
very large vacuum contribution in the initial entangled 
state, it is evident that this state is not highly entangled. 
However, by implementing the IPS scheme, the vacuum 
term is basically split in two new eigenstates (|01) and 
1 10)) leading to a large increase in the entanglement. By 
subtracting two photons (2PS scheme) rather than a sin- 
gle photon, the balancing of the eigenstates is destroyed 
for low photon numbers, and thus the entanglement is 
lower. However, it is clear that by implementing the 
displacement-enhanced scheme (D2PS), the weights be- 
tween the eigenstates are re-balanced, thereby producing 
strong entanglement. 

We now investigate the performance of the protocol un- 
der dissipation. This is of interest as in most protocols 
the entangled state is distributed in a network connected 
by lossy channels, thereby rendering the entangled states 
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FIG. 4: The logarithmic negativity is plotted as a function 
of the loss of the two channel. The attenuation is identical 
in the two channels. See caption of Fig. [2] for abbreviations. 
Inset: Purity and logarithmic negativity as a function of the 
displacement amplitude for the case of single mode squeezed 
vacuum (IMSV), a loss of 20% of the channels and a squeezing 
parameter of s = 0.3 (corresponding to an average photon 
number of 0.093). 
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in any practical realization impure. The logarithmic neg- 
ativity after distillation using the different strategies out- 
lined above for the single mode squeezed state is depicted 
in Fig. |4]for s\ = 0.3 (corresponding to (n) = 0.093) as a 
function of the channel attenuation. It is assumed that 
the two channels possess identical attenuations. As ex- 
pected, the LN decrease with the losses, but the effect of 
displacement-enhanced distillation is clear up to a loss of 
around 50%. We also find the LN as well as the purity 
as a function of the displacements as shown in the inset 

of Fig. m 

In the above study we have assumed that a single 
photon is perfectly removed from each of the entangled 
modes, and we have modelled this by the annihilation 
operator. This is however an idealization. In practice a 
photon is usually subtracted by reflecting a small part 
of the beam on an asymmetric beam splitter, and sub- 
sequently measuring the presence of a single photon in 
the reflected mode. The most practical detector for this 
purpose is an avalanche photodiode (APD) which is an 
on/off detector that discriminates between zero photons 
and some photons. Assuming that the reflectivity of the 
asymmetric beam splitter as well as the initial average 
photon number are low, the reflected beam will contain 
much less than a single photon on average, and thus the 
APD will effectively work as a single photon counter. On 
the other hand, for a finite refiectivity and a large initial 
squeezing, in some rare (although non-negligible) cases 
two photons will impinge onto the detector which cannot 



FIG. 5: Performance of the distillation protocol for realistic 
photon subtraction, a) Logarithmic negativity as a function 
of the initial average photon number, b) Success probability 
as a function of the average photon number. See caption of 
Fig. [2] for abbreviations. 



be discriminated from a single photon event. This causes 
an error. 

To simulate the realistic setup, we have used a more 
rigourous model in which the initial state (either ([T]) or 
([7])) is transformed through a beam splitter with a re- 
flectivity of 5%, displaced and finally measured with the 
projector 11 = 1 — |0)(0|. We use the package in Ref. [23] 
for the simulation. The results (both for the single mode 
and the two-mode squeezed state) are shown in Fig. [s] It 
is obvious from the figures that the trend of the LN as a 
function of the initial average photon numbers is identical 
to the trend for the ideal distillation scheme in Fig. [2^, 
but the amount of entanglement is lower in the former 
case. Another crucial parameter for characterizing the 
performance of the distillation protocols is the success 
probability. This is plotted in Fig. [5|d. At first it seems 
counter-intuitive that the success probability decreases 
when displacement is included, but it can be understood 
as follows. The rate of the initial photon detection in 
mode A is indeed increased when displacement is intro- 
duced. However, photon subtraction from a squeezed 
state increases its average photon number, so the dis- 
placement (which is experimentally implemented by ad- 
mixture of a coherent state) leads to a lower increase 
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in the photon number of mode B after the subtraction 
in mode A. Furthermore, the displaced photon subtrac- 
tion results in a state which has a small displacement in 
phase space (as opposed to the zero-mean of the initial 
state). The subsequent displacement of mode B before 
the photon detection there is opposite in direction of the 
state's displacement, leading to destructive interference 
in the detected mode. As a result of these two effects, 
the success probability of the second photon detection is 
considerably lower with displacement than without, out- 
weighing the increased probability of the first detection. 

In conclusion, we have theoretically investigated a 
displacement-enhanced distillation scheme of entangled 
states that are produced by a single squeezed mode. We 
have found that a simple Gaussian displacement opera- 
tion prior to photon subtraction increases the entangle- 
ment of the distilled state. Similar conclusion has been 
found for the two-mode squeezed state scheme, but in 
contrast to the previous proposals, the experimental re- 
alization of our scheme is much simpler as it does not 
require the control and phase locking of two independent 
squeezed beams. An experimental realization is therefore 
feasible with current technology 

While preparing this manuscript we were made aware 
of a related work [27] that builds upon our work in [22] . 
In Ref. [57], however, the robustness of the distillation 
scheme to a lossy channel was not investigated. 

We acknowledge financial support from the Danish 
Agency for Science Technology and Innovation. 
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